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ABSTRACT 

Pulsar braking torques due to magnetodipole radiation and to unipolar generator are considered, 
which results in braking index being less than 3 and could be employed to test the emission models. 
Improved equations to obtain pulsar braking index and magnetic field are presented if we deem that the 
rotation energy loss rate equals to the sum of the dipole radiation energy loss rate and that of relativistic 
particles powered by unipolar generator. The magnetic field calculated by conventional way could be 
good enough but only modified by a factor of ~ 0.6 at most. Both inner and outer gaps may coexist in 
the magnetosphere of the Vela pulsar. 

Subject headings: pulsars: general — radiation mechanisms: nonthermal 



1. INTRODUCTION 

Pulsar emission process is still poorly understood even 
more than 30 yr after the discovery. Nevertheless, it is the 
consensus of the researchers (e.g., Usov 2000) that primary 
pairs are produced and accelerated in regions (gaps) with 
strong electric field along magnetic line (-E|| ) and more 
secondary pairs (with multiplicity ~ 10^'^'*) are created 
outside the gaps (i?|| = 0), and that instability may be 
developed in the secondary relativistic plasma in or- 
der to give out coherent radio emission. Numerous models 
have been made concerning gap acceleration, whereas it 
is urgent to find effective way to test those specific and 
detailed models against observations. 

Because of observational difficulties, only braking in- 
dices n = U,Vl/Ci^ {Vl the angular velocity of rotation) 
of 5 young radio pulsars have been obtained observation- 
ally (Lyne & Graham-Smith 1998 and references therein, 
Camilo et al. 2000). They are: PSR B0531-f21 {n = 
2.51 ± 0.01), PSR B1509-58 {n = 2.837 ± 0.001), PSR 
B0540-69 (n = 2.2 ± 0.1), PSR B0833-45 {n = 1.4 ± 0.2), 
PSR J1119-6127 (n = 2.91 ± 0.05). Certainly, these ob- 
served indices include precious information on how pul- 
sars produce radiation. But all of them are remarkably 
smaller than the value n — ?> expected for the pure mag- 
netodipole radiation, according to which, the polar mag- 
netic field strength at pulsar surface, B, is conventionally 
determined by (e.g., Manchester & Taylor 1977), 



B 



sin a V Stt^R^ ' 



(1) 



where P — 2TT/il is the rotation period, / the moment of 
inertia, c the speed of light, R the pulsar radius, a the in- 
clination angle. B is singular (i.e., B oo) when a = 0°. 
Therefore, the i?-field derivation in this way is question- 
able and inconsistent since observation indicates n < 3, 
which means other processes do contribute to the braking 
torque. 

Indeed some efforts have appeared to find unusual 
torque mechanisms to understand the observed braking 
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index (see, e.g., Menou, Perna & Hernquist 2001, and ref- 
erences therein). An alternative effort, within the frame- 
work of "standard^ neutron stars and their magnetospheric 
emission models, is proposed in this Letter. We find that n 
and B derivation should generally depend on pulsar emis- 
sion models. Assuming that the orthogonal and aligned 
parts of magnetic moment are responsible to the dipole- 
radiation torque and to the unipolar-generator one, re- 
spectively, we obtain consistent equations for calculating 
braking index and magnetic field in the inner vacuum-gap, 
the space-charge-limited flow, and the outer gap models. 
We find that all of these models result in braking index 
n < 3, and in return the models could be tested for a par- 
ticular pulsar if its braking index and the inclination angle 
are observed. 



2. AN ASSUMPTION OF THE TOTAL ENERGY LOSS FOR 
ROTATION-POWERED PULSARS 

Pulsar broad-band emission depends essentially on a 
complete solution of the formidable well-defined magneto- 
sphere problem in relativistic electrodynamics and plasma 
physics, which, unfortunately, is still unknown hitherto 
(e.g., Mestel 2000). Nevertheless, the problem has been 
understood to some extent in two particular, i.e., the or- 
thogonal and aligned rotating, cases. 
Orthogonal Rotator An orthogonal rotator with mag- 
netic dipolar momentum fi± emits monochromatic elec- 
tromagnetic waves, the energy loss rate of which is = 
-^fi\n'^ ~ -6.2 X 10^7 BI2 Rl ^'^ ergs/s, where B12 = 
B/(10i2G), i?6 = i?/(106(cm)), and = BR^ /2. These 
low frequency waves are generally unable to propagate and 
should be absorbed in neutron star surroundings, and a 
larger amount of energy and the corresponding momentum 
could be pumped from neutron stars into their supernova 
remnants (Pacini 1967). 

Aligned Rotator The maximum potential drop in 
the open-field- line region by unipolar effect is (e.g., Ru- 

derman & Sutherland 1975) A$ = 



10« B12 Rl Vl^ cgse. 
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in charge depletion gaps, picking up energy in the gaps 
and angular momentum from the magnetic torque when 
streaming out. The angular momentum loss requirement 
(HoUoway 1977) can be satisfied if the charged particles 
can be "attached" to the magnetic field as far as near or 
out to the light cylinder. Two kinds of gaps are proposed 
to work in pulsar magnetospheres, which are called as in- 
ner and outer gaps. Various inner gaps are suggested, 
which depend on the binding energy of charged particles 
in pulsar surface, e.g., the vacuum gap model (Ruder- 
man & Sutherland 1975) with enough binding, the space- 
charge-limited flow model without any binding (Arons & 
Scharlemenn 1979, Harding & Muslimov 1998). The outer 
gap model was suggested to work near the null surface 
(e.g., Cheng, Ho & Ruderman 1986, Zhang & Cheng 1997) 
because the charged particles on each side of the surface 
should flow in opposite direction in order to close a global 
current in pulsar magnetosphere. It is thus obvious, as 
seen from above, that the energy loss is model-dependent 
for aligned rotators, which will be considered when calcu- 
lating pulsar braking indices and magnetic fields in the 
next section. Nevertheless, the energy loss rate of an 
aligned rotator, due to unipolar effect, could be written 
in the form of = — 27rrp • eg ■ A(/), if a gap has po- 
tential drop A(j> and the charge density in the gap is 
9 = CPgj ~ — 5.3CBi2fl cgse cm~^, where the po- 

lar cap radius Vp = R^/mfc ~ 5.77 x WrI'^Q}'^ cm. 
^ ~ 1 since g and Pgj are conventionally expected to be in 
a same order. 

The Assumption There are two schools of thought on 
the energy loss of an oblique magnetized rotator. One 
group opined that the magnetodipole radiation is the dom- 
inate mechanism of braking (e.g., Manchester & Taylor 
1977, Dai & Lu 1998, Lyubarsky & Kirk 2001), where no 
braking appears when a = 0. Another group suggested 
that pulsars' spindown dominates by a longitudinal cur- 
rent outflow due to unipolar generator (e.g., Beskin et al. 
1984), where Vl = constant if a = 90°. However, although 
there are two unseemly points when a = for the first 
school and when a = 90° for the second one, an interesting 
and strange thing, which is understandable in the next sec- 
tion, is that the derived physical parameters (e.g., B-field 
strength) are reasonable. We proposed that both energy 
loss mechanisms above, i.e., via dipole radiation and via 
unipolar generator, are expected to contribute the total 
braking torque of an oblique pulsar. Phenomenologically, 
for a pulsar with a total magnetic momentum fi = /2_l -|- /2|| 
(/xx — /isina, /^|| — /icosa), we could write the total en- 
ergy loss in the form of i? = cj^E'd + C||£'u, where c±_ and 
C|| are generally two functions of a indicating the contri- 
butions of those two energy loss mechanisms, respectively. 
Certainly c^_{a ~ 7r/2) = 1 and c\\{a = 0) = 1. An essen- 
tial and simple assumption to be employed in this paper 
is c_L = C|| = 1, since i? = _Ed + -Eu if and result 
independently in spindowns of E'd and -Eu, respectively. 
Therefore we have 



with 



(2) 



sin^ a 



3cos2aM 



A* 



~ sin^ a + 5.4 X IQ-'^R^^B^^ cos^ af^-^At/). 

3. BRAKING INDEX & ITS IMPLICATION 

The energy carried away by the dipole radiation (£^d) 
and the relativistic particles {E^ originates from the ro- 
tation kinetic energy, the loss rate of which is £^ = 10.(1. 
Energy conservation conduces towards 



2^ 
3^ 



(3) 



Based on Eq.(H), the braking index can be derived to be 
n = 6+ -. — = 



^df7' 



(4) 



which is not exactly 3 as long as rj is not a constant. If 
77 cx ri°, then n < 3 for a < (n > 3 for a > 0). For 
pulsars near death line, A0 ~ A^, i.e., the maximum 
potential drop A$ available acts on gap. In this case, 
r] — 1 + 2cos^a < 3, 77 = — 2 sin(2Q;)Q!. n < 3 if a gets 
smaller as pulsar evolves. For pulsars being away from 
death line, the potential drop A^ across an accelerator 
gap, which is model-dependent, is much smaller than A<I>. 
We discuss baking index in the following models, assuming 
that /2 (/i and a) and / are not changed for simplicity, since 
both observation (Bhattacharya et al. 1992) and theory 
(e.g., Xu & Busse 2001) imply that a pulsar's B-field does 
not decay significantly during the rotation-powered phase. 
The vacuum gap (VG) model The basic picture of 
vacuum gap formed above polar cap with enough binding 
energy was delineated explicitly in Ruderman & Suther- 
land (1975), where relativistic primary electrons emit 7- 
rays via curvature radiation in the gap. The gap potential 
difference A(/)^g = 4.1 x lO^p^^'^ B'^^^'^n^^'^ cgse, where 
the curvature radius'^ p ^ pe x 10^ cm. p ~ ^y^Rc/fl w 
2.3 X IO^Rq^^Q^^^^ for polar cap accelerators. We thus 
have A0VG ^ 9 2 x IO^^rI^'^ B^^^^'^fl-^^'^ cgse, r^^g ^ 

sin^ a + 4.96 x lO'^R'^^^'' B'^^"^ cos^ a^-^^/''. For vacuum 
gap where primary electrons emit 7-rays via resonant in- 
verse Compton scattering off the thermal photons (e.g., 
Zhang et al. 2000), the potential drop and the rj value 
are A^^G ^ i g ^ IQ^^ R^" B'^'''' n"/' cgse, 77^0 ^ 

sin2 a + 1.02 x lO^i^^^/^B-"'/^ cos^ a^-^^'^ 
The space-charge-limited flov^f (SCLF) model SCLF 
model works for pulsars with boundary condition of 
E\\ = at the pulsar surfaces. The previous SCLF (Arons 
& Scharlemann 1979) model has been improved to a new 
version (e.g., Harding & Muslimov 1998) with the inclu- 
sion of the frame-dragging effect. Though a simple and 
general analytical formula for all pulsar is not available in 
the Harding-Muslimov (1998) model, the potential drop 
could be well approximated in the extreme cases, regime 
I and II, which are defined as cases without or with field 
saturation^. In regime II case (i.e., the gap height be- 



ing larger than 



Zhang et al. 



(2000) obtained the 

1. But in this 



^Ruderman & Sutherland (1975) supposed there are multipole magnetic fields near pulsar surfaces, and they thus had pg 
paper we simply use dipole field lines for indication. 

*The definitions of regime I and II in Zhang & Harding (2000) have been misprinted (B. Zhang, 2001, personal communication) 
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potential drop, according to which rj values can be cal- 
culated. A</)SCLF ^ 7.1 x 10^ Rl^^n^/^ cgse, ryfCLF ~ 
sin^a + 38i?e:^/'*Sf2^cos2af7"^/4^ ^^le CR- 

induced SCLF models; A^^lic^s = 4.2 x 

^28/13^-9/13^18/13 ^g^^^ ^SCLF ^ ^-^2^ ^ 

2.3i?g"^"Bj"2^^/^^cos2ar2-**/i3, for the resonant ICS- 
induced SCLF models. In regime I, the stable ac- 
celeration scenario should be controlled by curva- 
ture radiation (Zhang & Harding 2000), Acpf^^^ = 
1.8 X 10"i?6/^Bi2'^'r!-i/^ cgse, rySCLF ^ ^^^2^ ^ 
9.8i?e""/^Bi'2'/' cos2 ar!-i5/^ 

The outer gap (OG) model For a self-sustaining 
outer gap, which is limited by the pair produced 
by collisions between high-energy photons from the gap 
and soft X-rays resulting from the surface heating by 
the backflowing primary pairs, the potential drop is 
A(j> = /^A$, where the fractional size of such outer 
gap / = 5.5Bi2'*/V26/2i (Zhang & Cheng 1997). / < 
1, which is satisfied for the five pulsars, if outer gap 
works. The rj value therefore can be calculated, A^'^^ — 
1.59 X 10i2i?3Bi2^/^17-io/2i cgse, r]'^^ ~ sin^ a + 8.6 x 
lO^Sis'^/^cos^ ar!-52/2i^ 

From these 77 values in different models, the braking in- 
dex can be obtained by Eq.(^. For typical pulsars with 
Re = 1 and B12 = 1, we compute the braking index n in 
each model, which is shown in Figjl]. It is obvious that 
n < 3 as long as inclination angle a < 90° in all of the 
models. Pulsars with small rotation periods tend to have 
n « 3. Also we can see from Fig. ^ or Eq.(^) that there is 
a minimum braking index 71(0; — 0°) for each model. In 
case of B12 = i?6 = 1, "-cr(" = 0°) = 0.86, = 
0°) = 1.14, n^'^ia = 0°) = 0.52, nf[;^l{a = 0°) = 1.25, 
nSCLF (a = 0°) = 2.38, nf^^^ia = 0°) = 0.86. 

We can not solve out magnetic field B by only Eq.(|3|) 
because rj = r]{a,^l). If a = 90° (or ?/ = 1), the solution 
of Eq.(||) results in Eq.(|]). In principal, Eq.(||) and ^ 
should be combined to find consistent B and a in case of 
braking index being known. However, because 1 < < 3, 
the magnetic field derived from Eq. (|l|) is good enough but 
only modified by a factor 1/y/rj £ (0.58, 1). 

Based on Eq.(||) and (^, the inclination angles of the 
five pulsars with observed braking indices are calculated in 
different models (see Table H). No solution of a is available 
for the Vela pulsar (PSR B0833-45) and PSR B0540-69 for 
the regime II SCLF(ICS) model since their braking indices 
are smaller than nff^^^a = 0°). This is consistent with 
the fact that these pulsars are young, and their gap heights 
are thus much smaller than Tp. 

Furthermore, we can determine whether a model works 
on a particular pulsar by comparing the calculated a in 
Table |l| with the observed a. Usually a can be derived 
by fitting the position angle curves of pulsars with high 
linear polarization in the rotating vector model (Lyne & 
Manchester 1988). For the five pulsars, only the inclina- 



tion angle of the Vela pulsar is obtained (~ 90°), however 
no a value in Table |l| tallies with this observation. There 
may be two possibilities to explain the discrepancy. (1). 
The braking torques due to the dipole radiation and to 
the unipolar generator should be treated and added in an 
other manner (e.g., Harding et al. 1999), rather than the 
way of ours. However, our treatment about the torques 
is reasonable, a further improvement of braking calcula- 
tion might not change substantially the results presented. 
(2). No model listed in Table |] can perfectly describe 
the actual accelerate situation of the Vela pulsar. The 
outer gap model explain well the high-energy emission of 
this pulsar, but could be still a partial description of the 
global magnetosphere. One possible picture is that both 
inner and outer gaps coexist in a pulsar's magnetosphere 
(Usov 2000), but the interaction between these two gaps 
and the pair plasma properties are still very uncertain. 
It is also possible that pair production process in strong 
magnetic and electric fields should be improved. For ex- 
ample, if B > O.lBc (-Be = 4.4 X lO^^G), 7-photons nearly 
along curved field lines convert into positroniums which 
could partially prevent the screening of (resultantly in- 
creasing the gap height and possibly having > 1), and 
therefore the energy loss increases significantly in po- 
lar cap models (Usov & Melrose 1996). Such an increase 
could result in a larger a in Table g (see Eq.(|)) since 
all magnetic fields of the five pulsar are very strong (near 
or greater than O.lSc)- In conclusion, further studies of 
testing emission models via braking index and of the theo- 
retical meaning of the test result would be interesting and 
necessary. 

4. CONCLUSION & DISCUSSION 

We have proposed in this Letter that the observed brak- 
ing index n < 3 could be understood if the braking torques 
due to the dipole radiation and to the unipolar generator 
are combined. The discrepancy between the observed in- 
clination angle and that derived from the six models of 
the Vela pulsar in Table |^ may call for improved pulsar 
emission models. In addition it is found that the magnetic 
field strength of a pulsar by conventional method could be 
a pretty good representation of the actual one. 

FigJ^ shows the variations of braking index n as func- 
tions of pulsar periods. Since pulsars spin down in their 
life, the curves in Fig.|l] represent the variations of n as 
functions of pulsar ages to some extent, n decreases as a 
pulsar evolves. However, the Johnston-Galloway's (1999) 
method to derive braking index can only be applied if n 
is constant during pulsar life. Therefore n can not been 
obtained by only P and P in principle. 
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Table 1 

The inclination angles (a) of the five pulsars derived from models 



Name (PSR) 


VG(CR) 


VG(ICS) 


OG 


SCLF(II,CR) 


SCLF(II,ICS) 
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B0531+21 
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2.9° 


5.0° 


2.1° 


1.6° 


33° 


B0540-69 
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6.3° 


5.2° 


1.8° 




36° 


B0833-45 


2.6° 


6.9° 


6.7° 


1.0° 
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11° 


8.4° 
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Fig. 1. — A set of calculated braJiing indices, as functions of rotation period, for six kinds of emission models. Pulsars are assumed to have 
polar magnetic field B = 10^^ G and radius R = 10® cm here. The inclination angles are chosen to be 0° (solid lines), 30° (dotted lines), 60° 
(dashed lines), and 90° (long-dashed lines). "CR" and "ICS" indicate curvature-radiation-induced and resonant iiivcrse-Compton-scattering- 
induced gaps, respectively. SCLF(Regime I): SCLF model without field saturation, SCLF(II,CR) and SCLF(11,1CS): SCLF model with filed 
saturation (Regime II). 



